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We explore both analytically and numerically the properties of doped t—J models on a class of 
highly frustrated lattices, such as the kagome and the pyrochlore lattice. Focussing on a particular 
sign of the hopping integral and antiferromagnetic exchange, we find a generic symmetry breaking 
instability towards a twofold degenerate ground state at a fractional filling below half filling. These 
states show modulated bond strengths and only break lattice symmetries. They can be seen as a 
generalization of the well-known valence bond solid states to fractional filling. 



I. INTRODUCTION 

Highly frustrated quantum magnets are fascinating 
and complex systems where the macroscopic ground state 
degeneracy at the classical level leads to many intriguing 
phenomena at the quantum level. The ground state prop- 
erties of spin 5 = 1/2 Heisenberg antiferromagnets on the 
kagome and the pyrochlore lattice remain still puzzling 
and controversial in many aspects. While the magnetic 
properties of the Heisenberg and extended models have 
indeed been studied for quite some time, the investiga- 
tion of highly frustrated magnets upon doping with mo- 
bile charge carriers has started recently 1-6 . Such interest 
has been motivated for example by the observation that 
in some strongly correlated materials, such as the spinel 
compound LiV204, itinerant charge carriers and frus- 
trated magnetic fluctuations interact strongly 7 . Further- 
more the possibility of creating optical kagome lattices in 
the context of cold atomic gases has been pointed out 8 , 
making it possible to " simulate" interacting fcrmionic or 
bosonic models in an artificial setting 9 . 

At this point we should stress that the behavior in a 
simple single-band model at weak and at strong corre- 
lations are not expected to be related in a trivial way. 
The weak coupling limit allows us to discuss the elec- 
tronic properties within the picture of itinerant electrons 
in momentum space based on the notions of a Fermi sur- 
face and Fermi surface instabilities (see e.g. Rcfs. 1,4). 
Considering for example the Fermi surfaces of a triangu- 
lar or a kagome lattice at half filling we do not find any 
obvious signature of the magnetic frustration present at 
large U. Although at weak coupling these systems do not 
seem to be particularly special, at intermediate to strong 
coupling the high density of low-energy fluctuations of 
the highly frustrated systems display characteristic fea- 
tures from which the physics of the frustrated system of 
localized degrees of freedom will emerge 2 ' 3 . 

In the following we study a class of highly frustrated 
lattices, the so called bisimplex lattices 10 , which are com- 
posed of corner-sharing simplices residing on a bipartite 
underlying lattice. We restrict ourselves to the triangle 
and the tetrahedron as the basic building blocks in the 



following. This class hosts lattices such as the kagome or 
the pyrochlore lattice. 

Our main result is the spontaneous symmetry break- 
ing between the simplices, that are located on different 
sublattices of the underlying bipartite lattice. This insta- 
bility takes place at the electron density n — 2/3 for the 
triangle based lattices like the kagome lattice and at the 
density n = 1/2 for the tetrahedron based lattices like 
the pyrochlore lattice and leads to a bond order wave 
(BOW) phase where the kinetic energy is staggered for 
the neighboring simplices. This instability is driven by a 
cooperative effect of the kinetic energy and the exchange 
interactions. 

The outline of the paper is the following: In Sec. II we 
introduce the lattices and the model and provide gen- 
eral arguments for the occurrence of the spontaneous 
symmetry breaking. These arguments are based on the 
analysis of the spectrum of isolated simplices and on 
a doped quantum dimcr model. We then report in 
Sec. Ill mean-field calculations for the kagome and the 
pyrochlore lattices which underline the symmetry break- 
ing tendency. In Sec. IV we present numerical results 
for the kagome lattice, that were obtained by exact diag- 
onalization (ED) and by the contractor renormalization 
algorithm (CORE). In Sec. V the ID analogue of the 
kagome lattice, the so-called kagome strip, is treated with 
two powerful methods that are available for ID systems: 
the density matrix renormalization group (DMRG) and a 
fcrmionic renormalization group (RG) and bosonization 
analysis. Finally we summarize and conclude in Sec. VI. 



II. GENERAL MOTIVATION 

The bisimplex lattices shown in Fig. 1 consist of corner- 
sharing simplices (triangles or tetrahedra) that are lo- 
cated on an underlying bipartite lattice. Therefore, we 
can separate the triangles and the tetrahedra into two 
different classes, which is visualized in Fig. 1 by a dif- 
ferent line-style (light and bold bonds). To refer to the 
simplices of a given class we call them "up" - and "down"- 
simplices, and we use the same terminology to distinguish 
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FIG. 1: (Color online) Four different bisimplex-lattices. The 
kagome lattice (a) and the pyrochlore lattice (b) together with 
their lower dimensional analogues, the kagome strip (c) and 
the checkerboard lattice (d). The two types of corner sharing 
units ("up" vs "down") are distinguished by the line width. 
They correspond to the bond order wave symmetry breaking 
pattern occurring at n = 2/3 on the triangle based lattices 
and at n — 1/2 on the tetrahedron based lattices. 

the triangles, the tetrahedra or the lattice bonds. 

In this section we provide two different arguments why 
for correlated electrons on such lattices at the fractional 
filling with two electrons per "up" -simplex (one electron 
per simplex) a spontaneous symmetry breaking can be 
expected, resulting in a phase where the "up"-simpliccs 
differ from the "down"-simplices. For the kagome and 
the pyrochlore lattice the inversion symmetry is broken 
in this phase, whereas for the lower-dimensional ana- 
logues, the kagome strip and the checkerboard lattice, 
translational symmetry is broken. The following two ar- 
guments apply generally to all bisimplex lattices shown 
in Fig. 1. They are simple and illustrate the basic un- 
derlying physics. In the remaining sections we will pro- 
vide detailed numerical and analytical evidence for sev- 
eral of these bisimplex lattices, which show that these 
arguments provide the correct picture. In the following 
we use (unless otherwise specified) the t—J model to de- 
scribe the correlated electrons on these lattices. The t—J 
Hamiltonian is given by 

<ii) ° 

+ j y^(Sj ■ Sj - - uiUj), 

(ij) 

where the restriction to the subspacc of singly occu- 
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TABLE I: Classification of the eigenstates of the t—J model 
on a triangle and on a tetrahedron. The degeneracy is given 
in the form r x (2S + 1), where r is the dimension of the 
irreducible representation of S3, respectively 54, and S is the 
total spin of the state. The asterix denotes the states retained 
in the CORE calculations for the kagome system, see text. 



pied sites is enforced by the projection operator V = 
n.i(l — n i-\ n ii)- The hopping amplitude t is always cho- 
sen to be positive. A negative sign of t will most likely 
induce ferromagnetic tendencies at the fillings we are 
considering 11 . 



A. The limit of decoupled simplices 

To get a basic understanding of the effect of doping 
in highly frustrated lattices we first consider the limit 
of decoupled simplices by turning the couplings within 
the "down" -simplices off. To connect this limit with the 
uniform lattice we use the parameter a (0 < a < 1) 
that tunes the coupling strength of the "down" -bonds as 
(at,aJ), while the coupling strength on the "up"-bonds 
is kept constant at (t, J). For a < 1 the inversion sym- 
metry between the "up"- and the "down" -simplices is 
explicitly broken by the Hamiltonian. The eigenvalues of 
Ht-j and their degeneracies are listed in Table I for a 
single triangle and a single tetrahedron. For t > and 
J > 0, there is a non-degenerate state with N e = 2 that 
has the lowest energy of all states and is separated from 
the remaining spectrum by a finite gap. This state has at 
the same time the lowest kinetic energy (— 2t or —At, re- 
spectively) of all states and gains the maximal exchange 
energy (—J) for two spins. This state is not frustrated 
anymore because it minimizes the kinetic and the ex- 
change energy at the same time. After having revealed 
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this particularly stable state with two electrons on a sin- 
gle unit, we know that for a = the system with two 
electrons per "up" -simplex has a non-degenerate ground 
state and a finite gap to all excitations. In the tight- 
binding model the gap decreases with increasing a and 
vanishes exactly at a = 1. For interacting electrons, 
however, the strong correlations that arc induced by the 
constraint in the kinetic term of the t— J model and the 
frustration of the exchange term, that arises from the 
strong frustration of the lattice, are accompanied by a 
tendency to occupy rather local states (in our case, they 
still keep a substantial part of their kinetic energy) and to 
pair up the electrons in nearest neighbor singlets. There- 
fore, it is a reasonable possibility that the system even 
for a = 1 has a finite gap and breaks the symmetry be- 
tween the "up" - and the "down"-simpliccs spontaneously 
in order to profit from the bipartite and non-frustrated 
structure of the underlying lattice. Such an instability 
has the character of a bond order wave - i.e., modulated 
expectation values of the bond energies - and yields an 
insulating state which breaks inversion or translational 
symmetry, respectively. 

B. Doped quantum dimer model 

In the previous sections we considered the case a = 0, 
where the Hamiltonian itself is not invariant under in- 
version symmetry. In order to get some insight into 
the mechanism of spontaneous symmetry breaking, it is 
desirable to treat "up"- and "down"-simplices on equal 
footing. In the following, we present a simple but il- 
lustrative model with a symmetric Hamiltonian that 
breaks spontaneously the symmetry between the "up"- 
and "down"-simplices in the ground state. 

A close inspection of the wavefunction of the lowest 
energy eigenstate of two electrons on either a triangle or 
a tetrahedron reveals that it consists of the equal ampli- 
tude superposition of all possible positions of the singlet 
formed by the two electrons: 

i<j 

where the normalization J\f = s/3 for the triangle and 
J\[ = V6 for the tetrahedron. This wavefunction moti- 
vates us to design a simple quantum dimer model which 
on each triangle prefers the exact wavefunction described 
above. Such a Hamiltonian reads for example for the 
kagome lattice: 

H QDM = -i£[|^)<A| + |A><A| + |A)^|+h.c.] 

A 

-^[|^)(V| + |V)(V| + |V)(^|+h.c.](3) 

v 

where the Hilbert space consists of all coverings of the 
kagome lattice with N c nearest-neighbor dimers and N c 



monomers, N c counting the number of unit cells. This 
corresponds to the situation at n = 2/3 in the t— J model. 
The interpretation is simple: the antifcrromagnetic ex- 
change term favors all the electrons to pair up into sin- 
glets, while the kinetic energy term delocalizes the sin- 
glets on a triangle. The quantum dimer model for the 
tetrahedron based lattices are defined by letting a single 
singlet resonate on a tetrahedron. This simple model al- 
lows us to find the exact ground state on these lattices. 
The ground state is twofold degenerate and each state 
is the direct product of equal amplitude resonances on 
the same type of triangles/tetrahedra, either all "up" or 
all "down" . In such a situation each resonating dimer 
can independently fully optimize its kinetic energy. The 
argument has much in common with the reasoning for 
the close packed dimer model on the pyrochlorc lattice 
discussed in Ref. 12. 

Although this model is only a cartoon version of the 
real electronic system, it illustrates how the tendency of 
the electrons to form nearest-neighbor singlets obstructs 
the motion of the singlets between corner-sharing sim- 
plices, but within a given simplex an individual singlet 
can hop without obstacles and optimize its kinetic energy. 
The bipartite nature of the underlying lattice allows for 
the localization of the singlets on simplices without in- 
terference and triggers in this way the spontaneous sym- 
metry breaking. 



III. MEAN-FIELD DISCUSSION 

In this section we present a mean-field calculation for 
the kagome lattice and the pyrochlorc lattice. The mean- 
field discussion is particularly valuable for the pyrochlorc 
lattice, as due to its higher dimension it is less affected by 
fluctuations and is presently not treatable with numerical 
methods. We can show that in the mean-field analysis 
the spontaneous inversion symmetry breaking, discussed 
in the previous sections, is also the natural and leading 
instability. 

We start with discussing the properties of the nearest- 
neighbor tight-binding model on the kagome and the py- 
rochlore lattice, given by 

H a = -flN -tJ^Yl Yl < 4+ W a m , - c r+ I /a n ,<7> ( 4 ) 

where a G {], |} is the spin index and the indices m,n 
run from zero to the dimension of the lattice, d. Further- 
more, r is an elementary lattice vector connecting unit 
cells and the vectors a , . . . ,a d point to the vertices of 
an elementary triangle (tetrahedron) in the kagome (py- 
rochlore) lattice, ao = 0. N is the total electron number 
operator and /i is the chemical potential. In the following 
we will use units where t = 1 and we will always choose 
fi = — 1 for the kagome and /i = — 2 for the pyrochlorc 
lattice, which corresponds to two electrons per unit cell. 
Hq can be diagonalizcd in reciprocal space and can be 
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written as 

= ^Ljkmff' ( 5 ) 

k,m,a" 

with 

-£ko = £ki = 6c, 6cm = d + 1 for m > 1. (6) 

For the kagome lattice we have 

£k = v 71 + 8 cos(fci/2) cos(fc 2 /2) cos([fci - k 2 }/2) (7) 

with fc m = k • a m and we refer to Ref. 13 for an explicit 
formula for 6c f° r the pyrochlorc lattice. 

The three bands of the kagome lattice consist of one 
flat band and two dispersing bands. The dispersing 
bands arc identical to the bands of a honeycomb lattice. 
They are shown together with the density of states per 
unit cell and spin in Fig. 2. Note, that around the points 

^ D(0 
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FIG. 2: The kagome bands and the density of states per unit 
cell and spin. The energy is measured in units of t. 

K and — K the dispersion shows a Dirac spectrum, i.e., 
the bands 6c0 an d 6ci touch at these points with linear 
dispersion. For the given chemical potential the Fermi 
surface reduces to points at K and — K and the density 
of states vanishes linearly with 6 i.e., we have D(£) oc |£| 
for small £. 

The four bands of the pyrochlore lattice consist of two 
flat bands and two dispersing bands. The dispersing 
bands are identical to the bands of a diamond lattice. 
They are shown together with the density of states per 
unit cell and spin in Fig. 3. Note, that 6c vanishes along 




rxWLT X -4 -2 024 



FIG. 3: The pyrochlore bands and the density of states per 
unit cell and spin. The energy is measured in units of t. 

the lines connecting X and W. The density of states also 



vanishes linearly at zero up to logarithmic corrections, 
i.e., we have D(£) oc |£| log(|£|) for small 6 

Systems with this form of the density of states at the 
Fermi level are neither band-insulators nor normal met- 
als, therefore, they are sometimes called semi-metals or 
zero-gap semiconductors. Although they have an even 
number of electrons per unit cell and no fractionally filled 
bands, they have no energy gap at the Fermi surface. 
Fermi surface instabilities are suppressed in this situa- 
tion. There is no Cooper-instability that leads to an 
obvious breakdown of perturbation theory for arbitrar- 
ily small attractive interactions, as the particle-particle 
polarization function involves the convergent integral 
/ d£Z?(£)/2|£| at zero temperature. For the half filled 
honeycomb lattice it has been shown that the Coulomb 
interactions lead to non-Fermi liquid behavior and that 
strong enough Coulomb interactions lead to antifcrro- 
magnetic order and to the opening of a charge gap 14-17 . 
The situation in the kagome and the pyrochlore lattice at 
the filling considered here is different. Because the lat- 
tices are not at half filling, it is not obvious that even ar- 
bitrarily large U would enforce a charge gap (Mott insula- 
tor) and an antifcrromagnctic order would be hampered 
by the frustrated topology of the lattice. However, if we 
consider the triangles (tetrahedra) as the fundamental 
units of our lattice we obtain the honeycomb (diamond) 
lattice and the properties of this underlying bipartite lat- 
tice will be reflected in the ground state and provide a 
way to circumvent the frustration effects. 

Our goal is to study the electron-electron interactions 
of the H t -j-v Hamiltonian given by 

At-J-V = VB a V + JS i ■ S o + Vn i n J > (8) 

(«) 

We will show that both the exchange and the repul- 
sion term favor the bond order wave instability. As 
the projection operator V is difficult to handle in ana- 
lytic calculations, the projection is often approximated 
by a purely statistical renormalization of the Hamilto- 
nian with Gutzwillcr factors 18 . We obtain a renormalizcd 
Hamiltonian without constraints given by 

Hr = 9t Ho+Y, J 9J S < ' S J + Vn * n J ■ ( 9 ) 

(ij) 

The renormalization is given by the Gutzwiller factors 
g t = 25/(1 + S) and gj = 4/(1 + 8) 2 and S is the hole 
doping measured from half filling. Note, that the nearest- 
neighbor repulsion is not renormalizcd by a statistical 
factor. 

In the following we will determine the critical J and V 
for spontaneous symmetry-breaking in this model within 
mean field theory. Superconductivity is a possible way of 
spontaneous symmetry breaking. As it is an instability 
in the particle-particle channel, the relation £k = £-k, 
which is ensured by inversion and time reversal symme- 
try, plays an essential role. Concerning the symmetry of 
the order parameter, we can restrict ourselves to singlet 
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pairing in the spin sector, because the nearest neighbor 
interaction is antiferromagnetic, and to s-wave pairing 
in the orbital sector, because in this way we obtain a 
nodeless, even gap-function. 

Another possibility of spontaneous symmetry breaking 
is an instability in the particle-hole channel. Such insta- 
bilities tend to occur if a nesting condition of the form 
£k = — £k+q is fulfilled. In general, this condition is not 
ensured by basic symmetries and therefore instabilities 
in the particle-hole channel arc much more special than 
superconducting instabilities. In our case, the relation 
£ko = — £ki can be considered as perfect nesting with 
q = 0. Therefore, the relevant question is which one of 
the two considered instabilities is dominant in our sys- 
tem. In order to answer this question we consider the 
following single-particle Hamiltonian, 

Atrial = H + Aphffph + AppiJpp (10) 

where we have introduced the two quadratic Hamiltoni- 
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Note that only the density of states enters these formulas 
because we are restricting ourselves to q = instabilities. 
The system spontaneously breaks inversion (U(l)) sym- 
metry, if the coefficient of A ph (A 2 p ) in Eq. (14) changes 
sign. If we assume that only one of the parameters V 
and J is nonzero, we obtain the following expressions for 
the critical values: 



J ph 
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(ii) 



H, 



The idea is to calculate the expectation value of H r for 
the ground-state of -fftriai and to choose the variational 
parameters A pp and A p h such that this expectation value 
is minimized. In terms of the operators that diagonalize 
H we can express the pairing operators as 



Aph = ^ixkltoalk^+h.c, 

^PP = e km7-kmT7kml + h.C. 

k,m 



(12) 



The numerical values for J c and V c are given in TA- 
BLE II. One can see, that the tendency for inversion 
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TABLE II: The parameters for the kagome (K) and the py- 
rochlore (P) lattice. The critical values are given in units 
of t. The coefficient of A p h (bond order wave) in Eq. (14) 
is negative for J > J ph (V = 0) or for repulsive V > L c ph 
(J = 0). The coefficient of A pp (superconductivity) is nega- 
tive for J > J c pp (V = 0) or for attractive V < L pp (J = 0). 



with the relations 
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(13) 



For small values of A p h and A pp we can expand the 
ground-state expectation value of H r in terms of A pp and 
A ph . Using the Wick theorem, we obtain up to higher 
order terms 
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where AE is the deviation from the ground-state expec- 
tation value with A p h = A pp = 0. A is the number of 
unit cells, t = tg t , J = Jgj, b is the number of bonds in 



symmetry breaking is much stronger than the tendency 
for superconductivity in both lattices and that both the 
antiferromagnetic J and the repulsive V support the in- 
version symmetry breaking. The integral J p h is large be- 
cause the factor £q — £ 2 takes its maximum at £ = 
whereas the factor (£ + /i) 2 in the integral 7 pp is much 
smaller for small values of £. In other words, supercon- 
ductivity has the handicap that the potential is propor- 
tional to the dispersion ek, therefore it is small at the 
Fermi surface and is only finite due to the finite value of 
/j,. The nearest-neighbor repulsion is harmful for Cooper 
(particle-particle) pairing, as can be seen from Table II. 
In the particle-hole instability, however, two particles 
tend to form a singlet on every second triangle (tetra- 
hedron) on the kagome (pyrochlore) lattice. In this way 
the singlet is still mobile and keeps — dt of its kinetic en- 
ergy and at the same time reduces the nearest-neighbor 
repulsion energy from 4"L/3 (3V/2) to V on every sec- 
ond triangle (tetrahedron). On the triangles (tetrahedra) 



without a singlet, the expectation value of the nearest- 
neighbor repulsion is however still 4V/3 (3V/2). In the 
limit where the kinetic energy is negligible (t <g; V, J) 
also other phases may appear. It is therefore important 
to emphasize that a finite kinetic energy is necessary to 
stabilize the bond order wave, because this phase arises 
due to the interplay between the kinetic and interaction 
energy. 

The limit of large V was recently discussed in the con- 
text of LiV2 04 by Yushankhai et al. in Ref. 19 for the 
pyrochlore lattice with n = 1/2. The possibility of inver- 
sion symmetry breaking was not considered in that study. 
But if V is of the order of t, the optimization of the ki- 
netic and the repulsion energy can lead to a compromise 
which breaks the inversion symmetry. 

In the bond order wave phase that we found in this sec- 
tion for the kagome and the pyrochlore lattice, the "up"- 
triangles (tetrahedra) have a higher expectation value of 
the kinetic energy than the "down" -triangles (tetrahe- 
dra). Furthermore a gap proportional to A p h opens at 
the Fermi surface. Therefore, the system made a transi- 
tion from a semi-metal to an insulator. This transition is 
similar to the Peierls metal-insulator transition, where a 
half filled system lowers spontaneously its crystal symme- 
try in order to open a gap at the Fermi surface. Phonons 
or the elasticity of the crystal play a crucial role in the 
Peierls transition. In our case, as we showed, the transi- 
tion can be driven by a purely electronic mechanism in 
an infinitely rigid lattice. In reality, the crystal structure 
will always relax and in this way additionally enhance 
the transition. 



IV. NUMERICAL RESULTS FOR THE 
KAGOME LATTICE 

In this section we compare the predictions obtained in 
the preceding sections to numerical results for the kagome 
lattice. We will first discuss some exact diagonalization 
(ED) results for the kagome lattice at n = 2/3, which 
confirm the BOW clearly and then we will study the de- 
pendence of the excitation gaps on the parameter a us- 
ing the contractor renormalization (CORE) method. In 
essence the numerical results corroborate the analytical 
predictions on the presence of a bond order wave insta- 
bility for the kagome lattice. 
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FIG. 4: (Color online) Correlation function of the kinetic en- 
ergy (Eq. 18) of a 21 sites kagome sample at n = 2/3 and 
J/t — 0.4. The black, empty bonds denote the same refer- 
ence bond, the red, full bonds negative and the blue, dashed 
bonds positive correlations. The line strength is proportional 
to the magnitude of the correlations. 

defined as: 

CW(*,J'),(M)] = (Kin(i,j)Kin(M)) 

-(Kin(i,j)) (Kin(M)>, (18) 

where 

Kin(i, j) =~Y, 4,„Cj, a + h.c, (19) 

a 

and and (k,l) denote two different nearest-neighbor 
bonds of the kagome lattice, that have no common site. 
This correlation function has been calculated for all dis- 
tances in the ground state of a finite kagome sample with 
21 sites, containing 7 holes, at J/t = 0.4. The result is 
plotted in Fig. 4. The reference bond uniquely belongs 
to a certain class of triangles ("up" -triangle in our case). 
Based on the theoretical picture one expects the corre- 
lation function to be positive for all bonds on the same 
type of triangles and negative on the others. This is in- 
deed what is seen in Fig. 4. We have also calculated the 
same quantity for J/t = 1 and J/t = 2 and the bond 
order wave correlations (not shown) were becoming even 
stronger for larger J/t. In this respect the ED calcu- 
lations confirm the qualitative picture developed above, 
that the homogeneous t— J model on the kagome lattice 
at n = 2/3 has an intrinsic instability towards a sponta- 
neous breaking of the inversion symmetry. 



A. Exact Diagonalization results 

The analytical arguments presented in Sec. II and 
Sec. Ill predict a bond order wave instability at filling 
n = 2/3. In finite, periodic systems this instability can be 
detected with a correlation function of the bond strength, 
cither of the kinetic term or the exchange term. Here we 
chose to work with the kinetic term, but the exchange 
term gives similar results. The correlation function is 



B. CORE results 

We know from Sec. II A, where we studied the limit 
of decoupled units, that the system has a finite gap in 
all sectors for a = 0. In this section we show how the 
gaps of the different sectors of excitations depend on the 
parameter a. 

The CORE algorithm 20-22 , which is based on the exis- 
tence of strong subunits (triangles) that arc only weakly 
coupled, provides a suitable method to perform such a 
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study. This is certainly true for small values of a where 
the system is naturally divided into weakly coupled sub- 
units. For a = 1 the CORE algorithm imposes a certain 
bias, as "up"- and "down" -triangles are not treated on 
equal footing, however, we have good evidence from ED 
and from the mean field analysis, that the system itself 
has a strong tendency to break inversion symmetry spon- 
taneously and therefore we can expect the CORE results 
to be reliable even for a = 1. In fact, we compare the 
CORE results to the ED results for the smaller clusters, 
and the good agreement shows that the finite excitation 
gaps at a = 1 are not an artefact of the method. 

The CORE method extends the range of tractable sizes 
of finite clusters, based on a careful selection of rele- 
vant low-energy degrees of freedom. In order to apply 
this algorithm, the lattice has to be divided into blocks; 
here, we naturally choose the "up" -triangles. A reduced 
Hilbcrt space is defined by retaining a certain number of 
low-lying states on each block. The choice of the states 
to keep depends also on the quantities to be obtained. 
While for a ground state calculation fewer states already 
provide good results, one has to retain usually more 
states to calculate the excited states. Here we choose 
to keep the 4 lowest states in the 3-electrons sector, the 
7 lowest states with 2 electrons and the 2 lowest states 
in the 1-electron sector. These states are denoted with 
an asterix in Table I. This choice leads to a reduction of 
the local triangle basis from 27 down to 13 states, thus 
allowing indeed to perform simulations on larger lattices 
than would be possible by conventional ED. 

Then, by computing the exact low-lying eigenstates of 
two coupled triangles, we calculate the effective interac- 
tions at interaction range two for each value of a and we 
neglect longer range terms. Comparison to ED data on 
the smaller clusters shows that this approach gives very 
good results. 

The basic excitation gaps of interest in the present 
problem are the spin gap, the single particle gap and 
the two particle gap. These arc defined as follows: 




0.5 1 0.5 1 0.5 1 
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a 
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A s=1 =E(N e ,l)-E(N e ,0), 



(20) 



FIG. 5: (Color online) Excitation gaps of the t—J model on 
the kagome lattice at J ft — 1 as a function of the parameter a, 
which denotes the ratio of the intertriangle to the intratrian- 
gle couplings. The gaps are obtained by the CORE method 
for different sample sizes (and geometries). For the smaller 
samples (18,vl; 18,v2; 21) ED data is shown for comparison 
at a = 1. 



is a strong dependence of the gap curves on the spe- 
cific sample. Note that there is no discrepancy between 
ED and CORE results. The second phenomenon can be 
understood from the discretization of the finite size Bril- 
louin zones : indeed, the measured gaps directly depend 
on the distance between the closest point in the Brillouin 
zone to the corner of the zone, the if-point, which is 
the point where the gap opens in the mean-field picture 
(Fig. 2). The 18, v2 and the 27 sites samples both contain 
this specific point and differ only slightly in the values of 
the gap. Thus, our numerical data support the claim of 
a finite gap for all a £ [0, 1]. The strong dependence is 
at the same time also a hint towards a sizable dispersion 
of the excitations in this system. 



Ai p = \[E{N e + 1, 1/2) + E(N e - 1, 1/2)] - E(N e , 0), 

(21) 

A 2p = 1 [E(N e + 2,0) + E(N e - 2, 0)] - E(N e , 0), (22) 

where E(N e , S z ) denotes the ground state energy in the 
sector with N e electrons and spin polarization S z . 

We have determined these gaps on kagome finite size 
samples at n = 2/3 and J/t = 1 containing 18 to 27 
sites. Two different versions of samples with 18 and 24 
sites have been treated (vl and v2). The results are dis- 
played in Fig. 5. There are two main observations: (1) 
the gaps do not close for any a £ [0, 1], giving additional 
evidence for the proposed symmetry breaking; (2) there 



V. DMRG AND RG RESULTS FOR THE 
KAGOME STRIP 

In this section we study with two different methods the 
kagome strip, shown in Fig. 1 c), as the ID analogue of 
the kagome lattice. This lattice has been introduced in 
Ref. 23, where it was shown to share some of the peculiar 
magnetic properties of the 2D kagome lattice. We report 
extensive density matrix renormalization group (DMRG) 
calculations 24 , both for the t—J and the Hubbard model. 
Furthermore, we apply the renormalization group (RG) 
and the bosonization techniques to the weak coupling 
Hubbard model. Our results show, that for both, the 
Hubbard and the t—J model, and irrespective of the cou- 
pling strength the BOW phase is realized in the kagome 
strip at n = 2/3. 



( a ) -E«T< c i,<T C j>) 




FIG. 6: (Color online) DMRG results for a L — 24 kagome ladder at J/t = 0.4 and n = 2/3. (a) Local bond strength deviation 
of the kinetic term. Red, full bonds are stronger (lower in energy) than the average kinetic energy per bond. Blue, dashed 
bonds are weaker than average bonds, (b) Local bond strength deviation of the exchange term. The color pattern are the same 
as in the upper panel. The thickness of the bond denotes the deviation from the average value per bond. Note that the pattern 
of the kinetic and the exchange term are in phase. 
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FIG. 7: (Color online) DMRG results for the alternation of 
the bond strength of the kinetic term and the spin exchange 
term as a function of inverse system size 1/L, for different 
values of J/t. 



A. DMRG results 

The kagome strip - being ID - offers the opportunity 
to perform DMRG simulations, thus allowing a rather 
detailed numerical study of large systems. We first dis- 
cuss the properties of the t— J model at n = 2/3 on this 
lattice, and then make a connection by investigating the 
Hubbard model at different values of U to the analyt- 
ical weak-coupling results of the following section. For 
both models we report numerical evidence for the pres- 
ence of the bond order wave instability for a large range 
of interaction strengths. 

In contrast to the periodic systems considered above 
within ED, the DMRG works most efficiently for open 
boundary conditions (OBC). In the present context this 
has the additional advantage that for even length L of 



the strip only one of the two degenerate ground states 
is favored, and we can directly measure the local bond 
strength. For the purpose of illustration we show the 
local bond strengths for a system of L = 24 in Fig. 6. 
The upper panel shows the difference of the local kinetic 
energy with respect to the average, while the lower panel 
shows the local expectation value of the spin exchange 
term, using the same convention. 

The calculated pattern resembles the schematic picture 
drawn in Fig 1 c). In order to address the behavior in 
the thermodynamic limit we measure the bond strength 
alternation, i.e., the difference between the expectation 
values of the operators Kin(i,j) and Si ■ Sj in the mid- 
dle of the system for different lengths L and values of 
J/t. The scaling of these quantities is shown in Fig. 7. 
The finite size corrections are rather small and all the or- 
der parameters extrapolate to finite values, irrespective 
of the value of J/t. Note that even for the case J/t = 
there is both a finite alternation of the kinetic energy 
and the magnetic exchange term. The alternation of the 
magnetic exchange energy is roughly the same for all val- 
ues of J/t. The alternation of the kinetic energy however 
is increased with growing J/t ratio. 

Next we address the question of the excitation gaps 
in the symmetry broken phase. The theoretical picture 
predicts an insulating state with a finite gap to all exci- 
tations above the two-fold degenerate ground state. We 
calculate the single particle charge gap and the spin gap 
defined in equations (21) and (20), respectively. The cal- 
culated gaps are shown in Fig. 8. The finite size gaps 
are extrapolated to L = oo with a simple quadratic fit. 
All gaps extrapolate to a finite value, in agreement with 
the predictions. The charge gap more or less follows the 
increase of the alternation of the kinetic energy shown 
in Fig. 7, i.e., the gap is roughly multiplied by a fac- 
tor three going from J/t = to 2. The behavior of the 
spin gap is mainly driven by the fact that it scales with 
J/t. Note that even in the case J/t = the spin gaps 
seem to remain finite. It will be an interesting question 
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FIG. 10: The kagome strip bands and the density of states 
per unit cell and spin. The energy is measured in units of t. 



B. Weak-coupling discussion 



FIG. 8: (Color online) DMRG results for the spin gap and 
single particle gap for J/t = 0, 0.4, 1, 2, as a function of inverse 
system size 1/L. 



We consider in this section the weak-coupling Hubbard 
model with the interaction 




nfinity 



FIG. 9: (Color online) DMRG results for the kinetic energy 
alternation for Hubbard kagome strips at n — 2/3 of length 
L = 32 and L = 48. The modulation is non-monotonous as a 
function of U/t and shows a maximum around U/t ~ 10 ~ 15. 



to characterize the precise nature of the charge and spin 
excitations. This will be left for a future study. 

The weak coupling RG calculations in the following 
section are performed for Hubbard onsite interactions. 
Although we expect the behavior of the t— J model and 
the Hubbard model at large U to be similar, we have 
explicitly calculated the alternation of the kinetic energy 
for the Hubbard model as a function of U /t. The re- 
sults displayed in Fig. 9 show that this quantity has a 
maximum around U/t ks 10 ~ 15, and interpolates be- 
tween the exponentially small order parameter at weak 
U/t and the result for the t— J model at J = 0, which 
corresponds to U = oo. These results therefore suggest 
that for the particular case of the kagome strip the weak- 
coupling phase is adiabatically connected to the strong 
coupling limit. 



Int 



(23) 



on the ID kagome strip. We map this weak local 
Coulomb repulsion on the original kagome strip on an 
effective interaction for the underlying half filled two- 
leg ladder. We will show in the following that for 
the one-dimensional kagome strip, applying the RG and 
bosonization techniques that were developed for the half 
filled two-leg ladder, we find in fact a charge density wave 
(CDW) instability in the effective model for the two-leg 
ladder which corresponds to a BOW instability on the 
kagome strip. 

The tight-binding bands of the kagome strip with 
fx = —t are shown in Fig. 10. The dispersing bands are 
the same as the bands of a two-leg ladder. The flat band 
originates from states that are trapped within one rhom- 
bus. The density of states has square-root singularities 
at ±i, ±3£ and a delta-peak at it. The Fermi surface is 
given by the 4 points ±fcpi and ±/cf2, where &fi = t/3 
and &F2 = 27r/3. There is a finite density of states at 
the Fermi surface. The kagome strip can be viewed as 
a kagome lattice tube, i.e., a kagome lattice with finite 
width and periodic boundary conditions. In order to see 
that the bands in Fig. 10 are in fact a cut through the 
kagome dispersion shown in Fig. 2, one has to shift one 
of the dispersing bands by n. This difference arises be- 
cause our notation is chosen to emphasize the similarities 
of the kagome strip to the two-leg ladder. 

In contrast to the kagome and the pyrochlore lattice, 
the density of states at the Fermi surface is finite for the 
kagome strip and wc therefore expect qualitative changes 
in this ID system even for weak interactions. We perform 
a weak-coupling RG and bosonization analysis for the 
kagome strip, and we show that the bond order wave 
instability is already present for arbitrary weak coupling. 
In this section we will only present the results of this 
analysis and refer to Appendix A 1 for further details. 

We derive an effective interaction for the two-leg lad- 
der, that corresponds in weak coupling to the local 
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Coulomb repulsion on the kagome strip (23). In this 
derivation we can drop terms that involve the high en- 
ergy states of the flat band and focus on the states in the 
dispersing bands. We denote the annihilation operator 
of these states by 7k, a = lk,i,a where k is the momen- 
tum along the strip and i € {1, 2} is the band index. If 
we rewrite the Hamiltonian Hi nt in terms of these new 
operators we obtain the interaction. 

ki...k 4 

where the prime over the sum restricts the sum to mo- 
mentum conserving k-values. For weak interaction we 
can replace kj in <7ki...k 4 by (fcpi,, it) and we obtain the 
simple expression 

flki...k4 = (6i li2 5i 3 i 4 +6i li3 Si 2ii + ^ 1 i 4 5i 2 i 3 )/6, (25) 

where q = k\ + ki — k% — k^. 

The effective interaction (24) can now be expressed in 
terms of left and right moving currents and in this way we 
find the initial values for the RG equations of the two-leg 
ladder. The integration of the RG equations with these 
initial values converges to an analytic solution that was 
identified by bosonization techniques as a charge density 
wave solution 25 . This means that the operator 

°cdw = \Y1 7fe,i, CT 7fc+7r,2,<7 + 7*+ w ,a 1 a7fc 1 i,«, ( 26 ) 

k,a 

acquires a finite value. The bond order wave order- 
parameter on the kagome strip is given by the expec- 
tation value of an operator Obow- 

The operators Ocdw and Obow transform identically 
under all symmetries of the system and, therefore, they 
describe the same phase. 

In addition, Ocdw is the effective operator on the two- 
leg ladder for Obow, i.e., if one does the same substi- 
tutions as we did for deriving Eq. (24) one sees that 
Obow - *■ Ocdw, if one chooses the right prefactor in 
the definition of Obow- 

We have shown that the bond order wave instability 
that is expected to occur at rather strong interactions 
according to the arguments of the preceding sections, is 
in fact already present in weak coupling for the one- 
dimensional kagome strip. Together with the DMRG 
results of the preceding section we provide convincing 
evidence, that the BOW phase on the kagome strip is 
present both in the Hubbard model for all values of U 
and in the t— J model. 



VI. CONCLUSIONS 

In summary we have studied the occurrence of a bond 
order wave instability in several bisimplex lattices. We 



provided evidence that this instability occurs quite gen- 
erally in these lattices at the fractional filling of one elec- 
tron per simplex (two electrons per "up" -simplex), if the 
correlations - i.e., antiferromagnetic nearest-neighbor ex- 
change and/or nearest-neighbor repulsion - are strong 
enough. 

In weak-coupling the physical properties of the system 
are dominated by the dimensionality of the lattice, by its 
fermiology and by the density of states at the Fermi en- 
ergy. We show that in the intermediate coupling regime, 
where the kinetic and the interaction energies are compa- 
rable, at the filling with two electrons per "up" -simplex, 
the physical properties of these highly frustrated lattices 
are dominated by local states on the simplex. The bi- 
partite and corner-sharing arrangement of the simpliccs 
allows the creation of isolated or only weakly interact- 
ing simplices with low energy by spontaneously break- 
ing inversion or translational symmetry. This knowledge 
provides a good starting point for series expansions or 
further CORE calculations. 

The magnetic interaction and the chosen sign of the 
dispersion leads to a tendency to form nearest-neighbor 
singlets and nearest-neighbor repulsion leads to a ten- 
dency to avoid configurations with more than two elec- 
trons per simplex. If the underlying lattice is bipartite 
the system finds a way to satisfy both tendencies simul- 
taneously by localizing singlets on every second simplex. 
This localization leads only to a partial loss of the kinetic 
energy, because the singlets can still delocalize within the 
simplex. It is the cooperation between the kinetic and the 
interaction energy which stabilizes the bond order wave 
state. Note, that the bond order wave instability does 
not lead to an inhomogeneous charge distribution on the 
lattice. 

The bond order wave states, which we find on the 
different lattices, provide a natural generalization of 
the well-known valence bond solid states (e.g. dimerized 
phases, plaquette phases) found in many frustrated spin 
models to situations away from half filling where a de- 
scription in terms of spin variables only breaks down. 
The density is still a rational fraction, but n = 2/3 in 
the kagome and kagome strip case while n = 1/2 in the 
pyrochlore and the checkerboard case. Approximately 
these states are direct products of singlets on triangles 
or tetrahedra, similar to the conventional picture of a 
dimerized phase. In contrast to the phases at n = 1 how- 
ever, the present instability involves a cooperative effect 
of both magnetic exchange and kinetic energy. 

An interesting task is to study the properties of a 
lightly doped bond order wave phase. It can be assumed 
that the bond order wave order parameter decreases 
rather quickly with doping. However, it is conceivable 
that away from the commensurate filling the bond order 
wave order parameter coexists with a small supercon- 
ducting order parameter. This phase would at the same 
time break lattice symmetries and the U{\) gauge sym- 
metry and would be therefore similar to a supersolid. A 
closer investigation of these issues has however be left for 
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further studies. 

In general, we conclude that the bond order wave in- 
stability occurs quite generally in bisimplex lattices and 
for physically reasonable models and interaction param- 
eters. We show that doping frustrated spin models can 
lead to new phases. Our study may also be a step towards 
the understanding of the interplay of frustrated mag- 
netic fluctuations and itinerant charge carriers, which 
could play a role for example in the unconventional heavy 
Fcrmion material L1V2O4 7 or in Na x Co02- 
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APPENDIX A: RG ANALYSIS 

1. kagome strip 

The tight-binding Hamiltonian for the kagome strip 
with periodic boundary conditions is given by 



^ = -*£££[< 



(Al) 



r— 1 v— ±1 o 



~t~ C r,0,<j{ C rM,a + C r+v,v,o) + nj 



where a g {Til} are the spin indices and N is the number 
operator. The chemical potential, /i, will be fixed to —t = 
— 1 in the following. It is convenient to introduce Fourier 
transformed operators 



ice {±1,0}, (A2) 



where the fc-sum runs over the L k- values in [— tt, tt). We 
can write this Hamiltonian in a diagonal form 

3 

Ho = £ £ £ Cfc.j ll,i, a l k ,l,vi 

k 1=1 cr 

and obtain the energies 

£ fe ,i = 1-2 cos(fc), £ M = -1 - 2 cos(fc), £ M = 3. (A3) 
and the operators 





1 


lk.2,a 


\/2/3 fe 


7k,3,a _ 





Pk -p k 

a k \f2 a k 
1 -V2a k 1 



Ck-l 
Ck,0 
Cfe. + l 



(A4) 



with a k = \/2cos(fc/2) and (3 k = \J~T~- i-af and k € 

[-7T,7T). 

The local Coulomb interaction introduced in Eq. (23) 
can be written as 



H 



4ee 



Int 



„t 



t 

£ L^i * 2 C k 1 ,x,-\^k2,x,l L -k 3 ,x,l^k4,x,-\- 

fcl.../S4 2:— — 1 

(A5) 

The sum over the momenta ki . . . ki is restricted, such 
that q = ki+k2—k 3 —k 4 is a multiple of 2tt. Note, that the 
appearing phase factor is important to determine the sign 
of the Umklapp scattering processes correctly. We obtain 
the effective low-energy Hamiltonian (24) from Eq. (A5) 
by doing the substitutions 



Cfc,±,, 



Cfc.O,. 



7fe,2,, 



(A6) 



These substitutions rules are obtained from Eq. (A4) if 
we set k = fcpi in the first row and k = kp 2 in the 
second row and drop the third row in the matrix of the 
transformation. 

In this way we can map the weak-coupling Hubbard 
model on the kagome strip on an effective weak-coupling 
model on the two-leg ladder. The problem of a weak- 
coupling two-leg ladder has been extensively studied by 
renormalization-group and bosonization techniques 25-27 . 
We will adopt here the notation of Lin et al. in Rcf. 25. A 
general weak interaction can be conveniently expressed in 
terms of left and right moving currents. Dropping purely 
chiral terms the momentum-conserving four fcrmion in- 
teractions can be written as 



+ f 1,-1 Hi, -1 1, 1 j — fij^Rii ' Jiji- 



(A7) 



To avoid double counting we set fa = 0. Furthermore, 
the symmetry relations /12 = /21 (parity), bi 2 = 621 (her- 
micity), and bn = 622 (only at half filling) hold. We have 
therefore six independent coefficients. For our interaction 
we find the values 

4% 1 = bZ 1 = U, 4b p 12 = bl 2 = j, 4/f 2 = /f 2 =|. 

In addition we have Umklapp terms given by 



n ~ a ij 1 Rij 1 Li] "n 1 



f T t 



12 = U 21 an( i U 



h.C. 



P - n,° 



with = 0, 
Here we have the values 



and 



0, u p 12 



U_ 

12' 



'11 



24' 



Integrating the RG equations with these initial values 
shows that the solution converges to the analytic solution 
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of the RG equation where all coupling constants except This solution was identified by bosonization techniques 
for b\ x and b p lx diverge with fixed ratios given by as a charge density wave solution (CDW) solution. 

/f 2 = -\ll2 = K 2 = -\bl 2 (A8) 



-uf a = -2< 2 = -2u p n =g>0. 
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